Low Bias Negative Differential Resistance in Graphene Nanoribbon Superlattices 
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We theoretically investigate negative differential resistance (NDR) for ballistic transport in semi- 
conducting armchair graphene nanoribbon (aGNR) superlattices (5 to 20 barriers) at low bias volt- 
ages VsD < 500 mV. We combine the graphene Dirac Hamiltonian with the Landauer-Biittiker 
formalism to calculate the current Jsd through the system. We find three distinct transport regimes 
in which NDR occurs: (i) a "classical" regime for wide layers, through which the transport across 
band gaps is strongly suppressed, leading to alternating regions of nearly unity and zero transmission 
probabilities as a function of Vsd due to crossing of band gaps from different layers; (ii) a quantum 
regime dominated by superlattice miniband conduction, with current suppression arising from the 
misalignment of miniband states with increasing Vsd; and (iii) a Wannier-Stark ladder regime with 
current peaks occurring at the crossings of Wannier-Stark rungs from distinct ladders. We observe 
NDR at voltage biases as low as 10 mV with a high current density, making the aGNR superlattices 
attractive for device applications. 

PACS numbers: 72.80.Vp, 73.22.Pr, 73.21. Cd, 68.65. Cd 



I. INTRODUCTION 

Graphene'^' has attracted much attention due to the 
possibility of new devices that may surpass their semi- 
conductor counterparts in both speed and reduced power 
consumption. This is expected due to the unique prop- 
erties of graphene, e.g., the high mobility of carriers, 
which can lead to high current densities, and the tun- 
ability of the bandgap. Additionally, building devices 
on the surface could facilitate optical absorption and 
emission. Particularly, negative differential resistance 
(NDR) is essential for many applications. ^ In semi- 
conductor resonant tunneling diodes ^ and superlattice 
structures, NDR is based on Fabry-Perot-type inter- 
ferences arising from the impedance mismatch between 
the various layers. These semiconductor NDR systems 
can also show interesting phenomena, such as intrinsic 
bistability due to charge accumulation. ' ' Pursuing the 
recent interest in graphene superlattices transport and 
thermal properties,^ '""^ it is a natural question to ask 
whether a graphene superlattice could exhibit similar fea- 
tures. 

The occurrence of Klein tunneling in graphene should 
be an obstacle to the NDR effect, as it gives a mono- 
tonically increasing contribution to the current. Narrow 
graphene nanoribbons overcome this limitation as the 
lateral confinement quantizes the Dirac cone into few- 
eV-wide bands. Tight-binding calculations show that it 
is possible to find NDR in these narrow nanoribbons at 
high bias voltages, 1-2 V.-^"'"' However, for integrated 
circuits a low bias mV regime is desirable to reduce power 
consumption. Low bias NDR can also be achieved in 
other graphene and bilayer graphene systems. '"^ 

In this work we consider an A^-barrier superlattice po- 
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FIG. 1: (Color online) (a) Metal-aGNR junctions and the 
modulated chemical shift Aep of the Dirac point across the 
aGNR (Refs. 2")-28). Aq (shaded regions) denotes the bar- 
rier and valley bandgaps. Here we consider square potentials, 
solid line. The dashed line shows the numerical results of 
Ref. ' . (b) Additional electrodes modulate the Dirac cone 
shift into a superlattice potential. The bias voltage Vsd is 
also shown, (c) Doped layers of a semiconductor superlat- 
tice can also modulate the local potential, (d) Schematic of 
the e — Vsd diagram of the source-drain transmission coeffi- 
cient showing crossings of the bandgaps Aq (black lines). The 
shaded regions delimit the energy range between the source 
fis ~ Aep and drain ^d = Ms — eVsD chemical potentials. 



tential on a semiconducting armchair graphene nanorib- 
bon (aGNR); Fig. 1. The electronic structure of the 
aGNR is a quantized Dirac cone, due to the quantiza- 
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tion of the transversal momentum A:„ , and can be metal- 
lic, = 0, or semiconducting, fc„o 7^ 0, depending 
on the width W of the nanoribbon; fc„p is the closest 
to zero transverse momenta. We choose W = 346ao, 
such that the aGNR is semiconducting with a bandgap 
Ao = 28 mcV; ao = 0.142 nm is the C-C distance. We 
use the transfer-matrix formalism to calculate the source- 
drain transmission coefficient Tsd across the superlattice 
potential along the aGNR, considering a finite bias volt- 
age VsD, revealing the electronic structure of the system; 
Fig. 2. The potential drop from source to drain follows a 
piecewise constant profile layer by layer; Fig. 1(b). The 
current is calculated within the usual Landauer-Biittiker 
formalism. 

We find low bias NDR at zero and room temperatures 
within three distinct physical regimes, (i) For wide lay- 
ers, the transmission across the bandgaps Aq is strongly 
suppressed, and nearly unity for energies away from the 
bandgaps. With increasing voltage, both barrier and 
valley bandgaps split and cross as shown schematically 
in Fig. 1(d), showing, at the coincidence region, a pat- 
tern of diamond-shaped structures with alternating re- 
gions of finite and suppressed transmission, thus leading 
to NDR. For narrow barriers resonant tunneling across 
layers become relevant, (ii) At zero bias, hybridization 
of resonant modes leads to minibands with finite, nearly 
unity, transmission; Fig. 2(b)-2(e). At very low voltages 
eVsu ^ 10 mcV (of the order of the miniband energy 
width) the resonant states misalign, thus breaking the 
minibands into off-resonance Wannier-Stark ladders with 
suppressed transmission. This gives rise to a single cur- 
rent spike near eVsD 10 meV. (iii) With increasing 
eVsD, rungs of ladders from distinct minibands cross and 
hybridize, showing a new set of resonant spikes in Tsd, 
Fig. 2(a), thus leading to current spikes and NDR. 



II. PROPOSED SYSTEM & MODEL 

The modulation of the Dirac cone into a superlattice 
potential can be achieved by different setups. It was 
shown that local charge-transfer effects between graphene 
and some metals (e.g., Al, Cu, Ag, Au, Pt) rigidly shifts 
the Dirac cone;^'^^""'^ Fig. 1(a). A series of metallic stripes 
over graphene can create the proposed superlattice po- 
tential; Fig. 1(b). Equivalently, the same structure can 
be obtained by selectively doping graphene regions in an 
alternate fashion. Additionally, the aGNR could be ar- 
ranged along the doped/non-doped layers of a cleaved 
semiconductor heterostructure; Fig. 1(c). Narrow sys- 
tems (< 400 nm) are desirable to keep transport ballistic 
at room temperatures. 

We consider low-energy excitations of graphene 
within the envelope function approximation, ■ i.e., the 
graphene Dirac Hamiltonian. The finite size of the 
nanoribbon requires vanishing wave functions at the 
edges, where for aGNR both A and B sublattices of 
the honeycomb lattice are present. This leads to van- 
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FIG. 2: (Color online) (a) Energy-voltage diagram of Tsd for 
TV = 5 barriers showing the evolution of the N — 1 hybridized 
modes [panels (b)-(d)] into Wannier-Stark ladders. Labels A, 
B, and C show the zero-bias hybridized modes in panels (a) 
and (d). Crossings of ladders' rungs from distinct minibands 
increases Tsd near Vsd = 30 and 50 mV. (b) Schematic of the 
modulated Dirac point (dashed line) , bandgaps Ao ~ 28 meV 
(gray area), and confined mode B' . In the transmission coeffi- 
cient Tsd across two barriers (o = & = 50 nm) (c), the confined 
mode B' shows up as a resonant spike near 230 meV. For (d) 
TV = 5, and (e) TV = 20 barriers the confined modes hybridize 
into TV — 1 spikes, building up a miniband. Similar resonances 
lead to minibands at energies away from the bandgap region 
Ao. 

ishing boundary conditions for the envelope functions at 
these edges. The validity of these boundary conditions 
is discussed in Ref. . Within this description, the elec- 
tronic structure of an aGNR is a quantized Dirac cone, 
£ = shvf^y fc^ -f fc^ . Here s = ±1 for the conduction 
and valence bands, V{ 10^^ nm/s is the Fermi veloc- 
ity, kx is the momentum in the longitudinal direction x, 
kn — WK /W — 47r/3ao is the quantized transverse mo- 
mentum with integer n, and W — 346ao ^ 50 nm. The 
fundamental gap is given by Aq — 2/iWf|fc„g| — 28 meV, 
with kna ^ -0.021 nm"^ 

To calculate the transmission Tsd = Tsd(£, fcn, Vsd) 
we use the transfer-matrix formalism, which relates the 
coefficients of the incoming and outgoing plane waves at 
the source and drain leads across the superlattice lay- 
ers (see the Appendix for details). We consider a piece- 
wise constant superlattice potential along the x direction. 
Figs. 1(b), through which the electronic structure of each 
layer is shifted by the local potential. In Figs. 2-4 we 
show Tsd only for fc„Q , as it contains the major contri- 
bution for the current in all investigated cases. 

The current density of Dirac electrons in graphene is 
given by j(r) = 4euf'0^(r)(T'0(r), where the factor of 4 
accounts for the valley and spin degeneracies, V'(r) is the 
envelope function spinor for the K or K' valley, and a = 
{ax,<Jy) are the Pauli matrices. Within the Landauer- 
Biittiker formalism, ' the current reads 

/oo 
Tsd(£, kn, Vsb) [/s(e) - /D(e)] de, (1) 
-00 

where /s(e) = {1 + exp[(£ - ^s)/fcBT]}"^ and /D(e) = 
/s(£ -I- Vsd) are the Fermi-Dirac distributions at the 
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source and drain, and fj,s is the source chemical potential. 
We truncate the sum over n to a few fc„ near fc„„ . 



III. RESULTS 

In Fig. 2(b) we consider a narrow graphene well with 
a = 50 nm and 5 — )■ oo. The solution of the graphene 
Dirac equation within the bandgap Aq region shows a 
confined state."'" This state corresponds to the resonant 
spike within the Ag region in Fig. 2(b) for two barriers. 
For N barriers the confined states hybridizes into — 1 
states, leading to minibands for large N; Figs. 2(c) and 
2(d). The minibands away from the Ag region occur due 
to reflections at each interface. For finite bias the mini- 
bands break into single resonant levels, Wannier-Stark 
ladders, as the confined modes from each layer misalign; 
Fig. 2(e). At the crossings of Wannier-Stark ladders from 
distinct minibands the transmission increases due to res- 
onant tunneling. 



NDR regimes 

To contrast distinct NDR regimes in our system, we 
discuss the current-voltage characteristics /-Vsd and the 
energy-voltage Tgo diagram for the following three cases. 
We compare five-barrier superlattices with (i) wide layers 
[Figs. 3(a) and 3(b)] and (ii) narrow layers [Figs. 3(c) and 
3(d)]. We then discuss (iii) a 20-barrier superlattice with 
narrow layers; Fig. 4. The dashed lines in the Tsd di- 
agrams delimit the zero-temperature window of integra- 
tion for /sD, defined between the source /zs = 230 meV 
and drain /^d — I^-s ~ VsB chemical potentials. 



1. "Classical" regime 

For wide layers, a = b — 100 nm, tunneling 
across bandgaps is strongly suppressed and the Tsd di- 
agram. Fig. 3(b), follows closely the diamond pattern in 
Fig. 1(d). For eVsD < Aep = 230 meV the current in- 
creases monotonically as the barriers bandgaps misalign. 
At the coincidence region, eVso ^ Aep = 230 meV, 
crossings of barrier and valley bandgaps lead to the dia- 
mond pattern of finite and suppressed Tsd- This alter- 
nation leads to the NDR near Vsd = 350 and 450 mV, 
in Fig. 3(a). The intensity of the NDR in this regime 
increases with the layers width, as the tunneling across 
bandgaps becomes more suppressed. The dashed curve 
in Fig. 3(a) is calculated with the limiting case where 
tunneling is completely suppressed across bandgaps, i.e., 
rpciassicai _ q j^cj-qss a bandgap, and 1 otherwise. Note 
the similarity of the dashed classical line with the exact 
/sD calculations in Fig. 3(a). 

For narrow layers, a = 6 = 50 nm in Figs. 3(c) 
and 3(d), the NDR due to classical regime is absent 
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FIG. 3: (Color online) Current and energy-voltage diagram of 
the transmission coefficient for five-barrier superlattice with 
a = b = 100 nm [(a) and (b)] and a = 6 = 50 nm [(c) 
and (d)]. The current- voltage characteristics are shown for 
T = 300 K and K. For wide barriers (a) and (b) the current 
follows closely the limiting "classical" case of TgD either 
across bandgaps, or 1 otherwise (dashed line). 



as it requires strong tunneling suppression. Interest- 
ingly, however, the Tsd diagram of a few narrow lay- 
ers clearly shows the evolution of the zero-bias mini- 
bands into Wannier-Stark ladders with increasing Vsd; 
Fig. 2(e). The Wannier-Stark ladders remain as individ- 
ual transmission spikes while there is an overlap of bar- 
riers (or valley) bandgaps. For eVsD > {N — l/2)Ao this 
condition is violated, and the tunneling across individ- 
ual bandgaps dominate. At the crossings of barrier and 
valley bandgaps, resonant effects are still visible in the 
Tsd diagram as stripes, corresponding to confined states 
between the overlapping band gaps; see Fig. 3(d) near 
e = -50 meV and Vsu = 400 mV. 



2. Mimband regime 

Considering a larger number of barriers, = 20 in 
Fig. 4, the aligned resonant modes hybridize into super- 
lattice minibands; Fig. 2. If /is is located within the 
miniband, at low biases the current is dominated by the 
transmission across these resonant modes. As the bias 
increases, the modes misalign breaking up the miniband 
into Wannier-Stark ladders. For five barriers. Fig. 2(a), 
the rungs of the ladders shows nonresonant transmis- 
sion peaks, and enhanced resonant transmission at cross- 
ings of the rungs (see Wannier-Stark ladder regime be- 
low). For 20 barriers, transmission through nonresonant 
rungs is strongly suppressed due to the larger number of 
bandgaps. At very low voltages. Fig. 4, the current ini- 
tially increases with Vsd as the transport occurs through 
the miniband. Near eVsD ^ 10 meV (of the order of the 
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miniband width) the miniband breaks up into the non- 
resonant rungs suppressing the current, thus resuhing in 
a pronounced current peak. 




Bias Voltage (mV) 

FIG. 4: (Color online) (a) Current-voltage characteristics 
and (b) Tsd diagram of a 20-barrier aGNR superlattice with 
a = 6 = 50 nm. In (a) the currents for and 300 K in the 
range < I^sd < 125 mV are multiplied by 6 and 10, re- 
spectively, for clarity. As the voltage increases the miniband 
near 230 meV, Fig. 2(e), breaks up as the resonant levels mis- 
align, leading to the pronounced spike near 10 mV for K. 
Near 50 mV the resonant levels return as resonant crossings 
of Wannier-Stark ladder rungs [see also Fig. 2(a)]. At the 
crossings Tsd increases, showing current spikes at both and 
300 K for VsD < 230 mV. For Vsd > 230 mV the current 
spikes arise from crossings of rungs at the coincidence region. 



3. Wannier-Stark ladder regime 

With increasing bias, rungs from Wannier-Stark lad- 
ders of distinct minibands cross. Fig. 2(a), creating 
new resonances through the superlattice layers. For 20 
barriers, where transmission from non-resonant rungs 
is strongly suppressed, the crossings show sharp Tsd 
stripes, e.g., at Vsd = 75, 110, 150, and 210 mV; 
Fig. 4(b). Each of these stripes, and others with lower 
contrast at smaller voltages, leads to current spikes in 
Fig. 4(a). The spikes broaden with increasing bias as 
the band gaps misalign. For eVsD > Aep = 230 meV, 
the crossings of broadened Wannier-Stark ladders from 
minibands near the barrier and valley bandgaps show 
diamond-shaped structures in the Tsd diagram, thus 
leading to a series of NDR spikes similar to the classi- 
cal regime. 



IV. CONCLUSIONS 

We have found that three distinct regimes can lead to 
NDR in semiconducting aGNR superlattices. (i) In the 
classical regime the NDR occurs as the bandgaps of dif- 
ferent layers cross with increasing Vsd- (h) For narrow 
layers and very low biases, eVso 10 meV, the trans- 
port is dominated by the resonant tunneling through the 
miniband, and the NDR occurs as the miniband breaks 
into Wannier-Stark ladders with increasing bias, (iii) For 
higher bias rungs of distinct ladders cross originating new 
resonances and current peaks. Interestingly, due to the 
high mobility of the carriers, we obtain low bias NDR 
peaks with high current densities. 

Final remarks 

The predicted NDR effects reported here are strictly 
valid for ballistic electronic transport through ideal 
aGNR superlattices. For relatively clean systems, how- 
ever, we expect detrimental effects such as those induced 
by disorder, impurities and structural defects"'^'^^'^^'''- to 
broaden the resonances in the I-V curves, thus possibly 
reducing the peak-to-valley current ratios. Interestingly, 
a recent calculation for the electronic transport through 
a single-barrier defined on a zigzag-terminated graphene 
nanoribbon shows evidence for a transport gap despite 
the gapless spectrum of the edge states of the system. ' ' 
Therefore, we expect that a superlattice defined on a 
zigzag graphene nanoribbon should exhibit transport fea- 
tures similar to those of the armchair case investigated 
here. The effects of edge irregularities, strong disorder, 
and interactions (even at the Hartree level) lie beyond 
the scope of the present work and deserve further study. 
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Appendix A: Transfer Matrix 

In this Appendix we detail the calculation of the trans- 
mission coefficient Tsd through the nanoribbon super- 
lattice via the transfer-matrix approach. We describe 
the potential across the system as piecewise constant; 
Fig. 1(b). In each layer the potential is a constant 
Vj = V^^^ - eVsnXj /L. The superlattice potential V^^^ 
is for valleys, and Vb = 230 mV for barriers (typical 
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value obtained from Refs. 25-28). The second term is 
the potential energy drop across the j^^ layers due to the 
electric field, where Xj is the coordinate of the center of 
the layer j, and L is the distance between the source and 
drain. 

The solution of the Dirac equation in each layer j (j = 
S and D for the source and drain, and an integer for the 
intermediate layers) is given by the plane- wave spinors ' 
ipj,n{x,y) = e*'^"^iy9j (x). For convenience we write the x 
component in a matrix form fj{x) = Mj{x)(j)j, where 
the components of the spinor (pj — {aj (3jY' denote the 
coefficients of the outgoing and incoming plane waves. 
The matrix Mj{x) is 



M,{x) 



ikl^'^x+iOj 



-ik^^ x — i9j 



(Al) 



The eigenenergies in each layer are £j.„ — ^ + 

Sjhv{\J (fci"'-')2 + fc2 ^ with Sj = +1 for the conduction 

band and sj = —1 for the valence band, ki?'' is the 
longitudinal momentum in layer j, fc„ is the quantized 
transversal momentum (conserved through the system), 

and 9j^n = tan~^(fc„/fci"''). 

The continuity of the spinors at the interfaces yields 
LPj{xjj+i) = ipj+i{xjj+i), where xjj+i is the position 
of the interface between the layers j and j + 1- Applying 
this matching throughout the system, we obtain a 2 x 2 
matrix equation connecting the coefficients from source 
and drain — Tm4>d, where Tu is the transfer matrix 



given by 



Tm ^Y[M-\x,,j+i)M,+,{x,,,+,). 



(A2) 



The definition of the reflected and transmitted waves 
depends on the sign of the electron energy at source sg 
and drain sd, such that the source and drain coefhcients 
are given by 



6S 



(1 r), if ss = 
(r 1), if ss = 

(t 0), if SD = 
(0 t), if SD = 



+1, 
-1, 

+1, 
-1. 



(A3) 
(A4) 



From the graphene Dirac Hamiltonian, the current 
density reads ji"*-* = 4:ev{ipj{x)a3;^j{x). At the station- 
ary regime the current flow at source and drain is the 
same, requiring the match j| = , from which we iden- 
tify the transmission coefficient Tgo, 



2 cos 9y) 
cos 9s 



(A5) 



This transmission coefhcient as a function of the energy 
reveals the electronic structure of the system, in which 
the confined modes in between the layers show up as 
resonant spikes and minibands; Fig. 2. 
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